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Abstract 

The distribution of return intervals of extreme events is studied in time 
series characterized by finite-term correlations with non-exponential decay. 
Precisely, it has been analyzed the statistics of the return intervals of extreme 
values of the resistance fiuctuations displayed by resistors with granular struc- 
ture in nonequilibrium stationary states. The resistance fluctuations are cal- 
culated by Monte Carlo simulations using a resistor network approach. It has 
been found that for highly disordered networks, when the auto-correlation 
function displays a non-exponential and non-power-law decay, the distribu- 
tion of return intervals of the extreme values is a stretched exponential, with 
exponent independent of the threshold. 

PACS numbers: 05.40.-a, 05.45.Tp, 02.50.-r 
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I. INTRODUCTION 



Fluctuations of prices in financial markets, wind speed data or daily precipitations in 
a given place for the same time windows are typically described by time series x{t) made 
by uncorrelated records^'^. In these cases, by considering the return intervals of extreme 
events associated with the overcoming of a given threshold q (i.e. the time intervals between 
two consecutive occurences of the condition x{t) > q), it can be seen that the are ex- 
ponentially distributed^'^'^. In other terms, the probability density function (PDF) of the 
distribution of the is given by^'^'^ : 

P,{r)^{l/R,)eM-r/R,) (1) 

where Rg is the mean return interval. Of course, the higher the threshold q, the bigger 
is the value of Rq. On the other hand, in the last years it has become clear that several 
other important examples of time series display long-term correlations^'^'^. This is the case 
of physiological data (heartbeats''^ and neuron spikes^), hydro-meteorological records (daily 
temperatures^'^'^°), geophysical or astrophysical data (occurrence of earthquakes^^'^^ or solar 
flares^^), internet traffic^ and stock market volatility^' records. Long-term correlated series 
are characterized by an auto-correlation function, Cx{s), decaying as a power-law: 

2 N-s 

Cx{s) =< XiXi+s >= Tt H ^i^i+s ~ (2) 

with exponent 7^ (correlation exponent) between and 1. In this case, the mean correlation 
time T, given by the integral over s of Cx{s) , diverges. The effect of long-terms correlations on 
the statistics of the has been first studied by Bunde et al.^ and the conclusions of this study 
were the following, i) The mean return interval Rg is left unchanged by the presence of long- 
terms correlations, ii) The distribution of return intervals becomes a stretched exponential: 

P,(r) = a exp[-(6 r/R^y] (3) 

where the two exponents 7 and 'jx were found to be equal, iii) The return intervals them- 
self are long-term correlated with a correlation exponent 7' ~ 7^;. As noted by Altmann 

2 



and Kantz in a very recent paper^, the result i), obtained by Bunde et al.^ by statistical 
arguments, can be identified with Kac's lemma^^ introduced in the context of dynamical sys- 
tems. The results ii) and iii) were obtained in Ref.^ on the ground of numerical calculations 
performed on long-term correlated and Gaussian time series generated by the algorithm 
described in Ref.^^. It must be noted that the result ii) only apphes to linear time series 
(i.e. to series whose properties are completely defined by the power spectrum and by the 
probability distribution, regardless of the Fourier phases) Apart from this restriction, 
the stretched exponential distribution of the seems to be a general feature in presence of 
long-term correlations in a time series^'^. It must be underlined that this fact has important 
consequences on the observation of extreme events: indeed it implies a strong enhancement 
of the probability of having return intervals well below and well above Rg, in comparison 
with the occurrence of extreme events in an uncorrelated time series. Furthermore, it must 
be noted that the stretched exponential distribution depends on the parameter 7 only, being 
a and b in Eq. (3) functions of 7, as shown in Ref.^. 

In this paper it will be analyzed the effect on the distribution of the Vg of the pres- 
ence of finite-term correlations with non-exponential decay, a situation which can occur 
in systems which are approaching criticality, where intermediate behaviors, consisting in a 
non-exponential and non-power- law decay of correlations can emerge"^. Precisely, in the fol- 
lowing section it will be studied the statistics of the return intervals of extreme values of the 
resistance fluctuations displayed by a resistor with granular structure in a nonequilibrium 
stationary state^^"^-*^. 

II. METHOD AND RESULTS 

The time series analyzed in this paper consist in the resistance fluctuations of a thin re- 
sistor with granular structure, in contact with a thermal bath at temperature To and biased 
by an external current /. The resistance values are calculated by using the stationary and 
biased resistor network (SBRN) modeP^"^^. This model describes a thin film with granular 
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structure as a two-dimensional resistor network in a stationary state determined by the com- 
petition between two stochastic processes, breaking and recovery of the elementary resistors. 
A broken elementary resistor (with resistance 10^ higher than the resistance corresponding 
to a normal elementary resistor) can be associated with a high resistivity region within the 
conducting material. Both processes are thermally activated and biased by the external 
current. The resistance of the network and its fluctuations are calculated by Monte Carlo 
simulations^^' Within this model, the level of intrinsic disorder in the network (average 
fraction of broken resistors in the vanishing current limit^^) is controlled by a characteristic 
parameter^": A = {Ed — ERj/ksTQ, where Ed and Er are the activation energies respec- 
tively of the breaking and recovery processes. The intrinsic disorder parameter A ranges 
between: A^m < A < Xmax, where Xmax corresponds to an homogeneous resistor (perfect 
network) and A = Xmm ~ to the maximum level of intrinsic disorder compatible with a 
stationary state of the network (stationary resistance fluctuations) 

In addition to this intrinsic disorder, the SBRN model considers also the presence of 
a disorder driven by the external current /. As a consequence, for a given value of A, 
and for a network of given size, nonequilibrium stationary states of the network exist only 
when I < Ib (breakdown threshold). For contrast, when I > Is the network undergoes 
an electrical breakdown associated with an irreversible divergence of its resistance^^"-*^^'^"^. 
For an arbitrary value of A this breakdown corresponds to a first order transition from a 
conducting to a non-conducting state of the network^^~^°. However, at decreasing A values, 
when A Xmm, the system becomes more and more close to its critical point^^^^'^. 

All the details about the SBRN model and its results can be found in Refs.^^~^^. However, 
it must be noted that this model provides a good description of many features associated 
with nonequilibrium stationary states and with the electrical instability of composites and 
granular materials^^^^'^^, including the electromigration damage of metallic lines^^. Finally, 
it must be underlined that, apart from the specific system described by the SBRN model, the 
method adopted here for generating the time series can be also viewed as a pure numerical 
algorithm for generating numerical series with different and tunable correlation properties. 
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FIGURES 
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FIG. 1. Normalized resistance fluctua- 
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FIG. 2. Auto-correlation function of the 



tions versus time (this last is expressed in time series in Fig. 1 (black-thick curve) . The 
simulation steps). The resistance values have grey, black-dashed and black-solid curves 
been normalized to provide a zero average represent the best-fit with a stretched expo- 



and a unit variance for x{t). 



nential, an exponential and a power-law, re- 
spectively. 



Then, by indicating with R the resistance of the network (expressed in Q) and with t 
the time (expressed in iteration steps), long R(t) time series (typically made of 1 -^ 2 x 10^ 
records) have been generated and analyzed for different values of A, of the external current 
and of the network size. Precisely, normalized series with zero average and unit variance 
have been considered: x{t) = {R{t)— < R >)/cr, where < i? > is the average value of the 
network resistance and a the root-mean-square deviation from the average. The analysis has 
been performed by calculating the auto-correlation function and the PDF of the x records, 
the return intervals Vg of the extreme values for different threshold q and their distribution 
Pq{r) (the values of q are expressed in units of a). For small A values (high level of intrinsic 
disorder), it has been found that displays a non-exponential and non-power-law decay. 
This behavior is different from that obtained for high A values (low level of intrinsic disorder) , 
where decays exponentially (consistently with the Lorentzian power spectrum reported 
in Refs.^^'^^'^°). By focusing on this situation, of interest for systems which are approaching 
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criticality, in the following of this section results will be shown concerning a network of size 
125 X 125, biased by a current I — Ib — 0.011 A and obtained by taking A — 0.33. 
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FIG. 3. Return intervals of extreme val- FIG. 4. Return intervals of extreme val- 
ues above the threshold g = 5 (in units of ues above the threshold g = 1.5 for the same 
root-mean square deviation) for the time se- time series of the previous figures. Here only 
ries of Fig. 1. the first 5000 intervals are shown. 

The x{t) time series is shown in Fig. 1 (only a small portion of the total number of 
records, = 1.6 x 10^, is reported in this figure). As evident from Fig. 1, the resistance 
fluctuations exhibit a strong non-Gaussianity and actually their PDF is well described by the 
Bramwell-Holdworth-Pinton distribution^^, as discussed in Refs.^^'^'^. The auto-correlation 
function of x{t) is reported in Fig. 2. The function significantly deviates from a single 
exponential and from a power-law while it is well fitted by a stretched exponential: 



C,{s) = A exp[(-s/6)T)] 



(4) 



with the following values of the fitting parameters: A — 1.23, b — 74.9 and 7 = 0.54. 
Many other functions have also been considered for the best-fit of the data. However, 
it has been found that the stretched exponential optimizes the best-fit procedure with the 
minimum numbers of fitting parameters. It must be remarked that a stretched exponential 
describes a behavior intermediate between a simple exponential decay (which is obtained for 
7 = 1) and a constant behavior (a limit of power- law) for 7^0. Moreover, the correlation 
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time corresponding to the expression (4) of is finite. 

The sequence of the A^^ return intervals of the values above the threshold g = 5 is plotted 
in Fig. 3 as a sequence of Nq impulses. The figure shows that a succession of very short 
return intervals (r^ <^ 100, the apparently empty portions of the horizontal axis) is followed 
by a succession of long intervals, indicating a strong clustering of the extreme events, a 
feature similar to that exhibited by the data of Bunde et al.® which instead concern long- 
term correlated records and very different from what observed for uncorrelated time series^^. 
Thus, the clustering of the extreme events is present even if the x records are not long-term 
correlated while they are characterized by a finite correlation time. This clustering of the 
extreme values persists also by lowering the threshold. This is shown in Fig. 4 which reports 
the sequence of the return intervals obtained for = 1.5. 
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FIG. 5. Double- logarithmic plot of the 



FIG. 6. The same as Fig. 5 but calcu- 



normalized probability density of the return lated after random shuffling the x-records. 
intervals for different thresholds q. 



Figure 5 displays the probability density function Pq{r) of the distribution of the return 
intervals as a function of Tq/Rq for different thresholds q ranging from —1.5 to 5.0 (the 
probability density has been normahzed to Pq{l) ). A double- logarithmic plot of Pq has 
been adopted because in this representation a stretched exponential function with exponent 
7 appears as a straight line of slope 7. Therefore Fig. 5 shows that the distribution of the 
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return intervals of extreme values of the x-series is well described by a stretched exponential 
and that the value of the exponent 7 is independent of the threshold g in a large range of 
g-values. This occurs even in absence of long-term corrrelations and in presence of a finite 
correlation time. This not-obvious result agrees with the conclusions of Altmann and Kantz 
formulated in their recent paper^. For comparison, Fig. 6 reports the probability density 
function Pq{r) of the distribution of the return intervals obtained after random shuffling the 
records of the same x-series: in this case 7 = 1, i.e. the distribution of the is exponential, 
as it must be for uncorrelated time series. 

III. CONCLUSIONS 

The distribution of return intervals of extreme events has been studied in time series 
with finite-term correlations. Precisely, it has been analyzed the distribution of return 
intervals of extreme values of the resistance fluctuations displayed by a resistor with granular 
structure in nonequilibrium stationary states. The resistance fluctuations were calculated 
by using the SBRN model based on a resistor network approach^'''"^^ . It has been found 
that for highly disordered networks, when the auto-correlation function displays a non- 
exponential and a non-power-law decay, the distribution of the is well described by a 
stretched exponential with exponent 7 largely independent of the threshold q. This result 
shows that the stretched exponential distribution describes the distribution of the return 
intervals of extreme events not only when long-term correlations are present in the time 
series^'^'^, but also when flnite-term correlations exist among the records, characterized by 
a non-exponential decay, a situation typical of systems which are approaching criticality. 
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